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ABSTRACT 

We investigate the range of covering factors (determined from the ratio of IR to 
UV/optical luminosity) seen in luminous type 1 quasars using a combination of data 
from the WISE, UKIDSS and SDSS surveys. Accretion disk (UV/optical) and ob- 
scuring dust (IR) luminosities are measured via the use of a simple three component 
SED model. We use these estimates to investigate the distribution of covering factors 
and its relationship to both accretion luminosity and IR SED shape. The distribution 
of covering factors (fc) is observed to be log-normal, with a bias-corrected mean of 
< log 10 fc >= —0.41 and standard deviation of a — 0.2. The fraction of IR luminos- 
ity emitted in the near-IR (1-5 (im) is found to be high (~ 40 per cent), and strongly 
dependant on covering factor. 
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1 INTRODUCTION 

There is clear evidence that all powerful AGN are sur- 
rounded by a geometrically thick distribution of optically 
thick material (Rowan-Robinson 1977; Lawrence & Elvis 
1982; Antonucci & Miller 1985; Edelson, Malkan & Rieke 
1987; Elvis et al. 1994; Richards et al. 2006). However de- 
spite nearly three decades of study the exact nature of the 
obscuring material has remained controversial. Recognising 
the difficulty in maintaining a smooth and geometrically 
thick rotating structure, Krolik & Begelman (1988) sug- 
gested that the material must be "clumpy" or filamentary 
in nature, a view supported by both recent resolution VLTI 
observations (Jaffe et al. 2004; Tristram et al. 2007) and the 
lack of strong 9.7 fim silicate emission in the mid-IR (Roche 
et al. 1991). 

This has motivated a large range of "clumpy" models, 
which assume a torus-shaped distribution of dust clouds sur- 
rounding the central AGN (Nenkova et al. 2002; Honig et al. 
2006; Nenkova et al. 2008, Schartmann et al. 2008; Honig 
et al. 2010; Stalevski et al. 2012). These models have been 
found to offer reasonable agreement with the observed prop- 
erties of large samples of AGN selected in a variety of ways 
(Alonso-Herrero et al. 2003; Ramos- Almeida et al. 2009; Mor 
et al. 2009; Landt et al. 2010; Alonso-Herrero et al. 2011; Deo 
et al. 2011) 
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However these models are almost exclusively phe- 
nomenological; they fail to prescribe a physical origin for 
the obscuring material, and are extremely difficult to pro- 
duce/maintain in real galaxies. This has long been seen as 
a weakness for "torus" models and has driven the devel- 
opment of physically motivated models, including; circum- 
nuclear star-bursts resulting from mergers (Wada & Nor- 
man (2002); Cattaneo et al. 2005; Kawakatu & Wada 2008), 
warped accretion disks (Greenhill et al. 2003; Lawrence & 
Elvis 2010; Hopkins et al. 2012) and accretion disk winds 
(Elvis et al. 2002; Elitzur & Shlosman 2006). 

While the coming generation of IR and mm-wavelength 
interferometers (e.g. MATISSE for VLTI, ALMA) should 
allow direct assessment of these physical models via direct 
imaging of the obscuring material, much can be gleaned from 
statistical studies of the unresolved properties of AGN. In 
particular the covering factor, i.e. the fraction of sight-lines 
to the AGN centre obscured by dust, is a potential dis- 
criminator of physical AGN models given sufficient statis- 
tics (Lawrence 1991; Treister et al. 2008; Lawrence & Elvis 
2010; Elitzur 2012; Hopkins et al. 2012). 

Until recently the limited availability of large-area mid- 
IR imaging has meant it was not possible to compile covering 
factors, traditionally probed via the relationship between re- 
processed dust and accretion disk emission, for large samples 
of AGN. The Wide-Field Infrared Survey Explorer ( WISE), 
has begun to revolutionise the study of AGN in the IR by 
performing a sensitive (~ 100 — 1000 times deeper than 
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IRAS) all-sky survey at near-to-mid IR wavelengths. Early 
work has already shown the potential of WISE to return 
very large samples of AGN; Mor & Trakhtenbrot (2011) pre- 
sented a WISE-hsseA study of 15, 928 quasars, a factor of 
~ 50 times the largest study possible with Spitzer (Richards 
et al. 2006). 

Here we make use of the WISE all-sky data release 
(Wright et al. 2010), in combination with near-IR obser- 
vations from UKIDSS (Lawrence et al. 2007), and optical 
photometry /spectroscopy from SDSS to study the distri- 
bution of covering factors, and its dependance on IR and 
UV/optical properties, of a statistical sample of the most 
luminous (Lbol > 10 46 ergs _1 ) quasars. A number of topics 
can be addressed with this large uniform dataset. In this pa- 
per our aim is to address only the simplest questions: what 
is the typical covering factor for luminous quasars? What is 
the distribution of these covering factors? do any existing 
models predict this correctly? and, finally, what is the range 
in mid-IR SED shapes seen? In ij2]we describe the construc- 
tion of our quasar sample, while ^describes our SED fitting 
approach. §3] presents our results on the distribution of cov- 
ering factors and the relationship between covering factor 
and IR SED shape. Finally i|6] presents our conclusions. 

Throughout we assume a ACDM cosmology with pa- 
rameters Qm = 0.3, Qa = 0.7 and Ho = 70kms _1 Mpc -1 . 



2 DATA 

The parent dataset for this study is the SDSS DR7 QSO 
catalogue (Schneider et al. 2010) as presented by Shen et 
al. (2011; henceforth Sll). This catalogue presents derived 
quantities, such as emission line fluxes, bolometric luminosi- 
ties and BH masses for 105,783 Type 1 quasars brighter than 
Mi — —22.0 with reliable redshifts. To avoid concerns about 
host galaxy contamination we restrict the luminosity range 
of our study to Lbol > 10 46 ergs -1 . In addition we only se- 
lect quasars with z < 1.5, as beyond this redshift we have 
little information about the rest-frame mid-IR (at z = 1.5 
the WISE 22 /mi band corresponds to rest-frame 8.8 /mi). 
This reduces our parent sample to 26,927 quasars. 

To build the WZSE-UKIDSS-SDSS (WUS) quasar sam- 
ple we cross-match these 26,927 quasars with first the 
UKIDSS Large Area Survey (LAS; Lawrence et al. 2007), 
then the WISE all-sky data release (Wright et al. 2011). 

The UKIDSS LAS is surveying 4000 sq. deg. of the 
sky in four near-IR bands; Y,J,H and K. The UKIDSS 
project is defined in Lawrence et al. (2007). UKIDSS uses the 
UKIRT Wide Field Camera (WFCAM; Casali et al, 2007). 
The photometric system is described in Hewett et al (2006) , 
and the calibration is described in Hodgkin et al. (2009). 
The pipeline processing and science archive are described 
in Irwin et al (2009, in prep) and Hambly et al (2008). We 
make use of the ninth data release (DR9) which has lim- 
iting magnitudes of 20.2, 19.6, 18.8 and 18.2 mags (Vega) 
at Y, J, H and K, respectively. The SDSS QSO sample is 
cross-matched to the LAS using a 2 arcsec search radius. Of 
the 26,927 quasars in our parent sample, 9230 have a coun- 
terpart in UKIDSS LAS with a detection in at least one 
near-IR band. 

The WISE final data release consists of imaging of the 
entire sky in 4 near-to-mid IR bands centred on 3.4, 4.6, 12. 



and 22 /mi to a depth of 0.04, 0.06, 0.5, and 3.2 mjy (3a). 
We consider WISE sources which are within 3 arcsec (one 
FWHM for WISE at 3.4 /mi) of UKIDSS-SDSS quasars to be 
reliable matches. This returns 9112 matches with a detection 
in at least one WISEb&nd; there remain 118 UKIDSS-SDSS 
quasars which do not have WISE counterparts. To ensure re- 
liable SED fitting in the IR we only select the 5281 quasars 
with reliable photometry in all four WISE bands (> 3er and 
CC_FLAGS equal to '0000'). For the 3831 quasars without 
four band WISE detections; 266 are excluded due to cor- 
rupted photometry (cc_FLAGS not equal to '0000') and are 
not considered in the following analysis; 3565 are genuinely 
undetected by WISE in at least one band (typically 22 /mi) 
and for these upper limits (3a) to the WISE photometry 
are estimated. These quasars, in addition to the 118 with no 
W'TS'iJcounterpars, will become important when we consider 
the selection biases of our study. 

In addition to the WUS quasars we also make use of 
SDSS quasars with IR photometry at 3.6, 4.5, 5.8, 8.0 and 
24 fim from Spitzer collated by Richards et al. (2006; hence- 
forth R06). We supplement the SDSS and Spitzer photom- 
etry presented by R06 with near-IR photometry at J and 
K from the UKIDSS Deep Extragalactic Survey (DXS) and 
12 /mi imaging from WISE. While the bulk of the R06 sam- 
ple is detected in the WISE bands we do not use the photom- 
etry in the 3.4, 4.6 and 22 /mi channels as it overlaps with 
the superior Spitzer photometry at these wavelengths. Of 
the 259 quasars present in the original R06 catalogue, 116 
have UKIDSS DXS and WISE 12 /mi counterparts within 
2 arcsec. While the R06 sample is much smaller than our 
new WUS sample, it does have the benefit of significantly 
deeper (~ 10 x ) mid-IR data from Spitzer. Again this will be 
important when considering the biases of the WUS sample. 



3 ESTIMATING AGN LUMINOSITIES 

Our aim is to use simple but reliable estimates of the accre- 
tion disk luminosity and the fraction of this re-processed by 
dust and emitted in the IR. We do this by fitting SED mod- 
els as described below. We do not require that the precise 
parameters of these model fits to be meaningful, as our SEDs 
have only a limited number of photometric bands. Rather, 
they are simply intended as an objective way to characterise 
the data. 

Sll present bolometric accretion disk luminosities for 
all of our WUS quasar sample, while R06 similarly present 
accretion disk luminosities (henceforth refered to as Lbol) 
for that sample, and we make use of these here. 

Fig. [T] shows Lbol as a function of redshift for both the 
WUS and R06 samples. The median luminosity and redshift 
of the WUS and R06 samples are comparable, with values 
of LZY S = 46.4 and z wus = 1.1 for WUS and medians 
of L^f = 46.2 and z Roa = 1.0 for the R06 sample. It can 
be seen that WISE is able to detect SDSS quasars above 
10 46 Lbol out to at least 2^1 and is not significantly biased 
to high Lboi when compared to the deeper R06 sample. 

For the IR luminosities we integrate the predicted flux 
density from the best-fit three component SED model. The 
three components of our SED model are; the accretion disk, 
a "hot" dust component, and an obscuring torus. No IR 
component from star formation is considered, meaning con- 
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Figure 1. Luminosity (-Lbol) vs. redshift for the WUS (red tri- 
angles) and R06 (blue filled squares) samples. The dashed black 
line shows the effective luminosity limit for quasars with a E94 
SED imposed by the i < 19.1 magnitude limit of the parent SDSS 
quasar sample. It is clear that the depth of the IR imaging used 
here is sufficient to probe the full range in z — -Lbol sampled by 
the SDSS QSO dataset. 



Figure 2. WISE colours [3.4-4.6] vs. [4.6-12] for WUS quasars 
(AB mags). Ovcrplotted are tracks taken by a quasar mid-IR SED 
(E94; solid line) and the local starburst galaxy M82 (dashed line). 
The colours of the WUS quasars show good agreement with the 
predicted E94 WISE colours, and are clearly separated from the 
predicted colours of M82. 



tamination in the mid-IR from star formation is effectively 
ignored. We justify this omission in two ways. Firstly, the 
observed WISE colours of our quasar sample are uniformly 
consistent with AGN dominated SEDs. In Fig. [2] we show 
the WISE colours [4.6-12] vs. [3.4-4.6] (AB mags) for our 
WUS sample, compared to the predicted location of quasar 
(E94) and the starburst (M82) SEDs, and the AGN colour 
selection "wedge" of Mateos et al. (2012). The WUS quasars 
show good agreement with the E94 quasar track, with the 
vast majority located within the AGN selection wedge. The 
WISE colours of WUS quasars are also clearly separated 
from the M82 track, very few of the WUS quasars lie close 
to the predicted colour evolution of M82 and in those cases 
the redshifts disagree significantly (i.e. z > 3 starbursts 
show similar colours to z ~ 0.1 quasars). The results of 
Fig. H strongly suggest that the mid-IR SEDs of WUS 
quasars are relatively unaffected by star formation. This is 
not surprising if we consider the luminosity of the WUS sam- 
ple. The minimum luminosity of our WUS QSO sample is 
ibol = 10 46 ergs -1 , for IR emission from star formation to 
be comparable to this would require a star formation rate 
(SFR) of ~ 450 Mq yr" 1 (assuming the translation between 
SFR and IR luminosity given by Kennicutt 1998). While 
this is feasible, indeed submm/mm observations of QSOs 
have revealed some host SFRs of 100-lOOOs M yr -1 (Isaak 
et al. 2002; Omont et al. 2003; Beelen et al. 2006; Hatzimi- 
naoglou et al. 2010; Dai et al. 2012), the typical SED for a 
star forming galaxy has an effective grey body temperature 
of ~ 40 K (compared to ~ 1500 K for AGN tori) and thus 
only a small fraction of the total IR luminosity is emitted at 
the mid-IR wavelengths probed by WISE. To quantify this, 
at the maximum redshift of our WUS sample [z < 1.5) the 
longest wavelength probed is ~ 9 fim. For the M82 SED used 
in Fig. [2] only 15 per cent of the IR luminosity is contained 
in the wavelength range 1-9 /im (compared to 60 per cent 
for the N08 AGN torus models considered here). 

At the minimum quasar luminosity we consider 
(10 46 Lboi), for star formation to contribute 0.2Lboi to the 
observable mid-IR luminosity in the WISE wavebands would 



require a SFR of 600 Mq yr 1 , while for the mean luminos- 
ity of the WUS sample ((Lboi) = 46.4), this would increase 
to ~ 2000 Mq yr" 1 . The question remains, are these SFRs 
common for quasar hosts at z < 1.5? Recently Spitzer and 
Herschel have provided far-IR observations for large samples 
of quasars for the first time (Shi et al. 2009; Hatziminaoglou 
et al. 2010; Bonfield et al. 2011; Dai et al. 2012). There is 
some evidence of a relationship between Lboi and SFR (as 
probed by Lfir), with the form of Lfir cx L^ 1 (Hatz- 
iminaoglou et al. 2010; Bonfield et al. 2011). Taking the 
results of Hatziminaoglou et al. (2010; Fig. 3) we can see 
that for submm detected type 1 AGN (which form only 15 
per cent of the total population) log 10 Lfir ~ 45 — 46 erg s - 
at Lboi = 10 46 ergs -1 , equivalent to a SFR of between 45- 
45OM0yr -1 . Given that these sources must represent the 
most star forming quasars (as ~ 85 per cent of quasars re- 
main undetected at Herschel wavelengths) the large SFRs 
needed to significantly affect the mid-IR SEDs considered 
here must be extremely rare in quasar hosts. 

We now describe each of the model components in de- 
tail. For the accretion disk component we assume a single 
template; the Elvis et al. (1994) radio-quiet quasar tem- 
plate from 0.05 < A < 0.7 /im extrapolated smoothly with 
a power-law of slope A/a cx A" 1 at A > 0.7pim (Elvis et al. 
1994; Hao et al. 2010). Emission from the Ha line at 6563A 
is also incorporated with a single Gaussian with a = 30A 
and a EW= 20. We also consider extinction of the accretion 
disk by dust, assuming SMC-like extinction curve (Prevot 
et al. 1984). In all cases the multi-component fit is made 
using the IDL routine MPFIT (Markwardt 2008). The IR- 
luminosities are calculated by integrating the best-fit SED, 
taking the hot-dust and "torus" components only in the 
range 1 < A < 1000 fim. 

For the hot dust we assume a simple blackbody de- 
scribed by two parameters; the amplitude and temperature. 
This component is added as previous work has found it dif- 
ficult to match the observed IR properties of quasars with a 
single clumpy torus (Mor et al. 2009) . To avoid degeneracies 
with the accretion disk and dust torus emission we restrict 
the temperature of the BLR hot dust to lie in the range 
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Table 1. Details of N08 model grid used to fit R06 quasars 



Parameter 


N 


Values 


i 


5 


0, 20, 40, 60, 80 


a 


5 


15, 30, 45, 60, 75 


Y 


4 


10, 30, 60, 100 


N 


1 


3, 6, 9, 12 


1 


3 


0, 1, 2 


TV 


5 


20, 40, 60, 80, 100 



100 < T < 2000 K. This upper limit is also consistent with 
the sublimation temperature for graphite grains (~ 1900 K; 
Barvainis 1987) 

Finally, the dust torus emission is characterised by one 
of three template SEDs selected from the clumpy torus mod- 
els of Nenkova et al. (2008; henceforth NO80 While the full 
range of the N08 models could have been used (N08 present 
over 100,000 different realisations of their model parame- 
ters) it is unlikely that the WISE 4-band IR photometry 
can distinguish between these models. Moreover, it is not 
computationally feasible to fit the full sample of over 5000 
WUS quasars with such a detailed model. Given that the 
intention here is to use the models to interpolate between 
the observed data points, not to interpret physical proper- 
ties of individual quasars, a small subset (i.e. three) of torus 
models should be sufficient. 

To confirm this, and to determine which template(s) to 
use, we selected a representative sample of 6000 N08 tem- 
plates and fit them to the smaller R06 quasar sample. The 
templates were selected on a representative grid in terms of 
the N08 model parameters; viewing angle (i), optical depth 
(ry), number of clouds (No), opening angle (a), power-law 
index of the radial distribution (q) and ratio of the exter- 
nal to internal radii (Y). The details of the N08 model grid 
are given in Table [1] The R06 sample is used as a testbed, 
rather than the WUS sample, as it has very similar charac- 
teristics (similar wavelength sampling, similar redshift and 
luminosity range) to the WUS sample, but is almost ten 
times deeper in the mid-IR, and so is not biased to mid-IR 
bright SED shapes. For each template on the grid we fitted 
the R06 sample using the full SED model described above. 

On the basis of these SED fits, we use a K-means clus- 
tering algorithm (Lloyd 1982) to group the R06 quasars by 
which N08 model provides the best-fit. The algorithm con- 
sists of iterating over two steps. In the first step each quasar 
is allocated to the group which offers the lowest \ 2 ■ In the 
second step the best model for each group is determined by 
finding the N08 template which offers the lowest mean x 2 - 
These steps are repeated until convergence (group assign- 
ments remain the same for two consecutive iterations). The 
initial templates for each group were chosen to be roughly 
representative of the spread in best-fit N08 parameter val- 
ues seen in the R06 samples. We run the algorithm assuming 

2, 3, 4 and 5 groups. The best single N08 torus model (i.e. 
one group) offers a mean \ 2 = 22.4. The mean % 2 across 
the R06 sample reduces from (x 2 ) = 18.7 in the two model 
case to (x 2 ) = 17.9, { X 2 ) = 17.6 and ( X 2 ) = 17.4 for the 

3, 4 and 5 model cases respectively. Given these results we 
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Table 2. Details of three N08 templates used to fit the quasar 
samples 
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justify the use of the three N08 templates to represent the 
WUS quasars; using more N08 models does not decrease the 
mean \ 2 considerably. 

The details of the three N08 models chosen by this pro- 
cess are given in Table [5] The parameter values we derive 
for the representative model set are noteworthy in two ways. 
Firstly, in all three templates our fitting calls for large values 
of Y, the ratio of outer to inner radius, and n the number 
of clouds, while requiring very low values of ry, the optical 
depth of the clouds. These values are at odds with other 
attempts to use the N08 models to describe the Spitzer IRS 
spectra of PG quasars (Mor et al. 2009), which typically 
call for lower values of Y and a and much larger values 
of ry. Secondly, while the first model would be seen as a 
Type 1 AGN the majority of the time, both the second 
and third model would rarely be seen as a Type 1 AGN; 
the accretion disk should be obscured ~ 80 and £ 90 per 
cent of the time given the parameter values used for the 
two models, respectively (Elitzur 2012). For these reasons 
it is unlikely that the parameter values from our fitting are 
physically meaningful, either because of our unsophisticated 
fitting methodology or the limited number of mid-IR wave- 
bands observed. However, our SED fits provide a meaning- 
ful estimate of the total IR luminosity and the gross shape 
of the IR SED by allowing us to interpolate between the 
mid-IR photometric points. Indeed, comparing the best fit 
SEDs from the N08 model grid to those derived just these 
three models the median \ 2 is found to increase (median 
X 2 = 11.1 for the grid vs. \ 2 = 12-8 for the three model 
case), but the derived quantities of interest (e.g. Lir and 
Lxspm) do not change significantly. To quantify this, for 
the IR luminosity we find the mean (fi) and standard devia- 
tion (a) of (Lf£ - Lfg 1 ) / Lgj to be fj, = 0.006 and a = 0.16, 
where Lf^ is the IR luminosity calculated from the best-fit 
to the full N08 grid, and Lf^ as the best-fit from the three 
templates. Similarly for Lis^m the mean (fj,) and standard 
deviation (a) of (Lf_ B/im - L? m 5M m)/ Lf_ 5)lm is found to be 
fi — 0.007 and a = 0.02. It should be noted that here (and 
throughout) Lir and Li_5 Mm do not include any contribu- 
tion from the accretion disk component. Assuming that the 
grid of N08 models completely describe the range of pos- 
sible "torus" shapes, these values represent the systematic 
and random error introduced by our adoption of a limited 
template set. 

For each WUS quasar we find the best fit combination 
of accrection disk, hot dust and torus components. We find 
the median X 2 of the WUS SED fits to be \ 2 = 19, slightly 
higher than the R06 sample despite the number of degrees 
of freedom being equal (Seven; thirteen bands less six free 
parameters). Of the 5281 WUS QSOs with reliable WISE 
photometry; 3600 are best fit by model (1) in Table [2] 1027 
by model (2) and 654 by model (3). Examples of five SED 
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Figure 3. Examples of SED fits to WISE-UKIDSS-SDSS (WUS) 
quasars. In each panel the solid line is our best fit SED, while the 
independent contributions from the accretion disk (dot-dashed 
line), broad line region hot dust (dashed line) and dust torus 
(dotted line) are also shown. The photometric data, de-redshifted, 
is shown as filled dots. Errors are lc. The panels show randomly 
chosen examples of (a) all with good fits (\ 2 < 20), (b) IR-faint, 
(c) IR-bright , (d) near-IR excess WUS quasars and (e) near- 
IR faint WUS quasars. The number of objects falling into each 
category is given in the top right of each panel. It can be seen 
that in all cases our model prescription does an adequate job of 
describing the observed SEDs of WUS quasars. 



fits which roughly span the range of SEDs seen in the WUS 
sample are shown in Fig. [3] 

It can be seen from Fig. [3] that, despite its simplicity, 
our SED modelling methodology has the ability to repro- 
duce the observed SEDs of a wide variety of WUS quasars 
with reasonable precision. Fig. [3] shows randomly chosen ex- 
amples of (a) all WUS quasars with good fits (\ 2 < 20), 
(b) IR-faint (covering factor less than one-fifth), (c) IR- 
bright (covering factor greater than one), (d) near-IR excess 
(Li_5 Mm /LiR > 0.55) and (e) near-IR faint WUS quasars 
(Li_5 Mm /LiR < 0.2). As well as illustrating the effective- 
ness of our SED fitting, these examples show that there is 
a significant range in the relative strength of the near-IR, 
mid-IR and accretion disk emission in the quasar popula- 
tion probed by WUS. 



4 RESULTS 

The ratio of the IR to the bolometric luminosity is com- 
monly interpreted as a "covering factor" , i.e. the fraction of 
sight-lines to the accretion disk which are obscured (Maolino 
et al. 2007; Hatziminaoglou et al. 2008; Rowan-Robinson et 



al. 2009). Fig. [4] shows the distribution of covering factors 
(LiB,/Lboi) for the WUS quasars. The mean (and standard 
deviation) of the observed covering factors is (log 10 /e bs ) = 
-0.34 and cr° O g S i0 ^ o = 0.19. To account for the signifi- 
cant fraction (~ 40 per cent) of quasars for which we only 
have limits to Lir, we assume a log-normal distribution for 
the data, and use the expectation-maximisation (EM) algo- 
rithm for maximum likelihood estimation of censored data 
presented by Wolynetz (1979). This results in a corrected 
mean and standard deviation of (log 10 fc) = —0.41 (i.e. 
{fc) = 0.3) and o"i ogl0 / c = 0.2. We can test our assump- 
tion of a log-normal distribution for the covering factor in 
two ways. First we can compare our bias corrected distri- 
bution to that seen in the smaller, but deeper, R06 quasar 
sample. The mean (and standard deviation) of the cover- 
ing factors in the R06 sample is (log 10 fc° 6 ) = —0.41 and 
c^g 6 f c = 0.19, almost identical to the corrected WUS val- 
ues. As second test we can use a Monte-Carlo simulation to 
predict what the observed distribution should be assuming 
the corrected one, and using the known limits of the WISE 
photometry. 

In each realisation samples are created with covering 
factors, redshifts and bolometric luminosities drawn from 
measured distributions. The detectability of each source is 
tested by assuming the mean mid-IR SED and predicting 
the observed WISE fluxes, assuming Gaussian errors con- 
sistent with the catalogued values. Sources are considered 
detected if they have WISE flux densities above the 3<r de- 
tection limits quoted in SJU subject to incompleteness levels 
consistent with the values presented in the WISE explana- 
tory materia^]. This process is repeated 100 times to produce 
statistically robust results. In our Monte-Carlo realisations 
we detect on average 61.9 per cent of sources, comparable 
to the actual detection rate of 59.6 per cent (5281/8846). 
The mean and standard deviation of the recovered cover- 
ing factors in the simulation is (log 10 fc m ) = —0.33 and 
"to™ f c ~ 0-19> m good agreement with the values for the 
observed distribution in WUS quasars. 

A potential variation is seen as a function of luminosity, 
with the mean log 10 fc decreasing by 0.11 across the two lu- 
minosity bins. Although this trend cannot be confirmed with 
the limited luminosity range of our sample, it is consistent 
with that seen in the analogous studies of WISE quasars 
by Mor & Trakhtenbrot (2011) and Calderone, Sbarrato & 
Ghisellini (2012). 

As our ability to constrain the shape of the IR SED 
with typically only 4-6 photometric points is limited, we test 
variations in the IR SED shape by comparing the near-IR 
(1-5 fxra) to the total IR luminosity (L 

1 — 5 fJ.m 

/Lir). Figure 

shows Li_5 Mm /LiR as a function of bolometric luminos- 
ity, IR luminosity and covering factor. While I/i-s^m/im 
appears relatively insensitive to Lboi and Lir, a strong cor- 
relation appears between Li_5 ^m/LiR and fc- To test the 
significance of this correlation we calculate the Spearman 
rank correlation coefficient (p), finding p — —0.5. Given the 
sample size this correlation has a << 1 per cent probability 
of occurring by chance. 
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Figure 5. Ratio of near-IR (1-5 (im) to the total IR luminosity (Li— 5 ^/Lu) as a function of; Lbol (left) , Lir (middle) and covering 
factor (right). Open triangles are all WUS quasars with reliable photometry, while filled triangles represent those with good SED fits 
(x 2 < 20). Filled symbols are colour-coded by the N08 model fitted (see Table [2}; with red representing model 1, green is model 2 and 
blue is model 3. Black contours enclose 90, 70 amd 50 per cent of the WUS population. A strong correlation between L\— 5 ^m/Lut and 
fc can be seen. 
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Figure 4. Distribution of covering factor (fc] defined as 
Lia/Lbol) f° r WUS quasars. The solid line shows the binned his- 
togram of fc's for WUS quasars with good (> 3<r) WISE detec- 
tions in all four bands. The dashed line shows the best-fit Log- 
Normal distribution to the full WUS sample, taking into account 
upper limits for those quasars without 4-band WISE detections. 
The top two panels show the distributions in bins of log 10 -Lbol; 
while the bottom two panels shows the distribution for all WUS 
quasars. Also shown in the third panel is the distribution for the 
R06 sample (blue dotted line). In the bottom panel the distri- 
bution of fc for Type 1 quasars is shown as a red long-dashed 
(dot-dashed) line for the warped accretion disk model of LE10 
uncorrected (corrected) for the selection limits of WISE. 



5 DISCUSSION 

5.1 The typical covering factor for high 
luminosity quasars 

The simplest result is that the mean observed value of 
log 10 fc =-0.41, i.e. the typical quasar is observed to have 
fc ~0.39. This is consistent with early SED studies such 
as those of Sanders et al. (1989) and Elvis et al. (1994), 
but seems to be inconsistent with the SDSS/Spitzer study 
of Treister et al. (2008), who found L 2 4/L hol ~ 0.04 at the 
highest luminosities. 

It is important to note that, under the assumption that 
the Type 1/2 divide is purely an orientation effect, our ob- 
served distribution of covering factors must be biased, such 
that pi(fc) = (1 — fc)p(fc), where pi(fc) is the probabil- 
ity of seeing a Type 1 AGN (quasar) with covering factor 
fc and similarly p(fc) is the probability of seeing any AGN 
with covering factor fc (Lawrence 1991, Lawrence & Elvis 
2010, Elitzur 2012). Interestingly, applying this correction 
gives us a median covering factor for all AGN of 0.64, in 
good agreement with the ratio of Type 2 to Type 1 AGN 
seen in the IR/radio selected samples (LE2010; based on De 
Grijp et al. 1992; Rush et al. 1993; Lacy et al. 2007). How- 
ever it is in serious disagreement with claims based on X-ray 
samples that at the highest luminosities the obscured frac- 
tion is only 5-10 per cent (Hasinger et al. 2008; Tueller et 
al. 2008), although this extra X-ray obscuration may orig- 
inate from dust-free gas clouds very close to the accretion 
disk (e.g. Risaliti et al. 2007; Elvis et al. 2012). 

5.2 The distribution of covering factors 

It is striking that there is a clear distribution of covering 
factors, but that the distribution is fairly narrow. The dis- 
persion in log 10 fc of 0.2 means that two-thirds of type 1 
quasars have fc in the range 0.25-0.61. In principle the 
shape of the covering factor distribution is a test of models 
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for AGN obscuration. In practice such models have not so 
far made such predictions. The exception is the simple tilted 
disk model of Lawrence & Elvis (2010; henceforth LE10), 
which predicts p(c) = | sin irfc (We ignore the other model 
proposed in LE10, the twisted disk, as it is known to give 
an unrealistic distribution of covering factors). 

We show this prediction for the tilted disk model in 
Fig. [4l taking into account the bias for Type 1 AGN, i.e. 
Pi(c) — |(1 — fc) sin nfc ■ This very simple model, with 
no adjustable parameters, comes close to the observed dis- 
tribution, but with predicted values slightly higher than the 
observed ones. This suggests that a detailed physical warped 
disk model may be able fit the observed distribution. It 
would be very interesting to see predictions from disk wind 
and star-burst disk models. 

To determine the effect of the WISE selection limits on 
the LE10 tilted disk model we performed a Monte-Carlo sim- 
ulation by assigning covering factors randomly selected from 
the LE10 tilted disk distribution to our real WUS quasars. 
This process was repeated 100 times and a prediction of the 
observed distribution made, excluding sources which would 
now be below the WISE detection limits and including a 
random error of 15 per cent on our measurements of the 
covering factor; this is shown in Fig. [4] Interestingly, our 
simulation returns roughly the right completeness, we would 
expect to detect 5185 quasars on average (similar to our ac- 
tual number of 5281), and a log-normal distribution, with 
the peak of the distribution offset by log 10 fc ~ 0.1 dex from 
the observed distribution. In order to bring our observed dis- 
tribution in line with these predictions our estimates of the 
covering factor must be systematically biased low by ~ 25 
per cent, or an additional absorber unassociated with the 
quasar must be present (e.g. galaxy scale dust; Elvis 2012). 

5.3 The relationship between covering factors and 
near-IR emission 

The results of Fig.[5]present us with two puzzles; is the frac- 
tion of IR luminosity emitted in the near-IR consistent with 
AGN models? and why is it a strong function of covering 
factor? 

Taking these questions in turn, the typical ratio of near- 
IR (1-5 fim) to total IR luminosity seen here is ~ 40 per cent. 
This is close to the maximum Li_s ^ m /£lR achievable within 
the N08 models, and is only possible with a limited range of 
model parameter values (e.g. Y < 30, N ^ 4, q ^ 1, r ^ 20). 
This discrepancy is not unique to our work, several previous 
studies attempting to explain the IR SEDs of AGN have 
found the need for extra "hot" dust components (Polleta 
et al. 2008; Mor, Netzer, & Elitzur 2009; Deo et al. 2011; 
Vignali et al. 2011). 

A potential reason for this excess may be the dust com- 
position. Thermal emission peaking at ~ 1 /im require tem- 
peratures which are greater than the sublimation tempera- 
ture for the silicate-type grains (T~ 1500 K) that make up 
the bulk of the dust assumed by the models (N08; Draine 
& Lee 1984). However some dust must survive very close to 
the accretion disk in order to explain the correlation be- 
tween the optical and near-IR variability seen in nearby 
AGN (Minezaki et al. 2004; Suganuma et al. 2006; Kishi- 
moto et al. 2007). While the single blackbody employed here 
is unlikely to be the correct model for the hot dust, the Wien 



side of the blackbody fit is well constrained by the UKIDSS 
and WISE near-IR data and hence we may interpret our fit- 
ted blackbody temperatures as a "maximum" temperature 
for the dust closest to the accretion disk. Encouragingly, the 
mean "maximum" temperature for WUS quasars is 1481, K, 
below the sublimation temperature for silicate grains, and 
in good agreement with the value (1400 K) found via a simi- 
lar analysis of PG quasars by Mor, Netzer & Elitzur (2009). 
Only 37 per cent of WUS quasars are found to need a hot 
dust in excess of 1500 K, suggesting that variations to the 
dust composition are not needed to explain the IR SEDs of 
most luminous quasars. 

The assumptions made by the "torus" model have a 
large effect on the near-IR luminosity. Models which assume 
a smooth, rather than clumpy, distribution of dust (e.g. Fritz 
et al. 2006) naturally result in a large fraction of the IR lumi- 
nosity emitted in the near-IR (see discussion in Vignali et al. 
2011). As pointed out by Vignali et al. (2011), the failure of 
clumpy models to reproduce the large near-IR luminosities 
seen here without additional hot components points to seri- 
ous problems with the underlying paradigm of the "clumpy 
torus". However smooth models are not without their own 
failings; both in terms of their stability (Krolik & Begelman 
1988) and their ability to reproduce the range of observed 
10 /im silicate emission/absorption (Nikutta, Elitzur & Lacy 
2009). 

The recent model of Stalevski et al. (2012) tries to in- 
clude aspects of both clumpy and smooth models by using 
a "two-phase" model; essentially a clumpy model embedded 
in lower-density, smoothly distributed, dust. Computing the 
near-IR to total IR luminosity for their template SEDqj we 
find that near-IR luminosity ratios similar to that seen here 
(e.g. Li_5 Mm /LiR > 0.4) are easily achievable for these mod- 
els. 

Finally, we consider the implications of the correlation 
between the ratio of near-IR to total IR luminosity and cov- 
ering factor. A simple solution is that the near-IR (i.e. "hot" ) 
and mid-IR (i.e. "warm") components of the SED come from 
physically distinct components. If the covering factor of the 
"hot" component was relatively constant across all quasars, 
while the "warm" component varied, it would produce a re- 
lation similar to that seen in Fig. [5] (as the "hot" and "warm" 
covering factors must add up to fc)- However, this would 
require an obscurer for the "hot" component physically dis- 
tinct from the "torus" which is needed to explain the mid-IR 
emission. Given the result shown in Fig. for this scenario 
to work the covering factor of the "hot" dust component 
must be ~ 5 — 15 per cent. Interestingly, this is close to the 
value assumed for the broad line region clouds (~ 10 per 
cent; Wyithe & Loeb 2002). 

Another explanation for the observed fc - Li_5 Mm /LrR 
relation, which does not require two distinct obscurers, is 
that it is a geometric effect. As the covering factor decreases, 
the maximum inclination at which a type 1 quasar would 
be seen increases. An increase in the inclination will mean 
direct sight lines to more of the "inner wall" of obscuring 
material closest to the accretion disk. 

Qualitatively it is possible to replicate this postulated 
relation between Li_5,j m /LiR and inclination in the N08 
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models, with a factor of ~ 2 increase in the ratio of near- 
IR to total IR luminosity from inclination 0° — 60° assum- 
ing model parameters describing a dense, compact torus 
(N £ 10, q k 2, Y < 10, a < 30). Given that the a and 
No parameters describe the covering factor (Elitzur et al. 
2012) this behaviour in the N08 models is consistent with 
our geometrically motivated solution; thin (low a), dense 
(high N, high q, low Y) tori should result in large near-IR 
to total IR luminosity ratios when seen at moderate incli- 
nations (;> 40°). Because the opening angle for the torus is 
fixed at 50° in the Stalevski et al. (2012) model a similar 
test for these models cannot be performed. 



6 CONCLUSIONS 

We have considered the optical to mid-IR properties of a 
sample of type 1 quasars selected from a combination of the 
WISE, UKIDSS and SDSS datasets. Using our simple SED 
modelling approach we estimate a number of quasar proper- 
ties; the IR- luminosity, the covering factor (from Lm/Lboi) 
and the IR SED shape characterised by the ratio of near-IR 
(1-5 /im) to total IR luminosity. Using these measurements 
we reach the following conclusions: 

(i) The distribution of covering factors (fc), defined as 
the ratio of IR to UV/optical luminosity, is found to obey 
a log-normal distribution. Once selection limits are taken 
into account the distribution is characterised by (log 10 fc) = 
—0.41 and o"i oglo f a = 0.2. These values agree well with other 
IR/radio estimates, but are well below that estimated from 
comparable X-ray studies. 

(ii) This distribution gives roughly the same shape as that 
expected for the tilted-disk model (Lawrence & Elvis 2010), 
although offset to lower covering factors by ~ 25 per cent 

(iii) A significant fraction (~ 40 per cent) of the total 
IR luminosity is emitted in the near-IR. It is difficult to 
replicate this behaviour with models in which all the torus 
material is in "clumps" (e.g. Nenkova et al. 2008). Smooth 
or two-phase models, which assume some clumpy material 
embedded in an otherwise smooth torus, show much higher 
near-IR to total-IR ratios, consistent with our observations. 

(iv) A strong correlation is observed between the ratio of 
near-IR to total IR luminosity and fc- This is interpreted 
as a geometric effect, as more of the hotter dust close to the 
accretion disk will be visible at the high inclinations possible 
in low fc quasars. 
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